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ABSTRACT 

The virialization shock around an X-ray cluster primarily heats the ions, since 
they carry most of the kinetic energy of the infalling gas. Subsequently, the ions share 
their thermal energy with the electrons through Coulomb collisions. We quantify the 
expected temperature difference between the electrons and ions as a function of radius 
and time, based on a spherical self-similar model for the accretion of gas by a cluster in 
an D = 1, h-so = 1 universe. Clusters with X-ray temperatures Tx = 4-10 x 10^ K, show 
noticeable differences between their electron and ion temperatures at radii ^ 2Mpc. 

High resolution spectroscopy with future X-ray satellites such as Astro E may be 
able to determine the ion temperature in intracluster gas from the width of its X-ray 
emission lines, and compare it to the electron temperature as inferred from the free-free 
emission spectrum. Any difference between these temperatures can be used to date the 
period of time that has passed since the infalling gas joined the cluster. 


Subject headings: galaxies: clusters of - cosmology: theory 


1. Introduction 

The standard technique for measuring the total mass in clusters of galaxies from X-ray 
observations assumes that the cluster gas is in hydrostatic and thermal equilibrium (see, e.g. 
Jones &: Forman 1992). Hydrodynamic simulations often maintain the assumption of thermal 
equilibrium, and show that while the gas is typically at hydrostatic equilibrium in the cores 
of clusters, its state is perturbed in their outer parts (Navarro, Frenk, & White 1995; Evrard, 
Metzler, & Navarro 1996). In this paper we address the question of whether the gas in these outer 
regions should actually be in thermal equilibrium as assumed. To see why the cluster gas might 
deviate from thermal equilibrium we should first examine how the gas obtains its high X-ray 
temperature. 

As fresh material joins the stationary cluster gas, it converts its infall kinetic energy into heat 
through a shock. Since most of the inertia of the infalling gas is carried by the ions, they acquire 
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most of the dissipated energy and heat to about twice the virial temperature of the cluster. The 
electrons, on the other hand, remain cold after the shock with their temperature slowly rising as a 
result of their Coulomb collisions with the hot ions. Thermal equilibrium will only be achieved if 
the temperature equilibration time between the ions and the electrons is shorter than the age of 
the cluster. In fact, the ratio between the equilibration time and the age of the Universe is of order 
unity, ~ (fIfe/i 5 o/ 0 . 06 )(T/ 10 ® K)^/^((5/178)“^, for the characteristic overdensity 6 and temperature 
T in cluster shocks, where fl;, is the baryonic density parameter and hso is the Hubble constant 
in units of 50 km s“^ Mpc“^. This implies that the electron and ion fluids might not achieve 
complete thermal equilibrium within the cluster lifetime. 

A difference between the electron and ion temperatures was invoked as a potential explanation 
for the increase in the inferred baryon fraction in the outer parts of A2163 (Markevitch et al. 
1996). The ASCA data show a sharp drop in the electron temperature beyond 2 / 159 ^ Mpc in this 
cluster. The gas distribution, determined by ROSAT, is well described by a /3-model. Assuming 
hydrostatic equilibrium, a dark matter distribution with the same profile as the gas was strongly 
ruled out by the temperature measurements, and models in which the dark matter fell off much 
more rapidly were only marginally allowed. If, however, the electron temperature is lower than the 
ion temperature, the actual gas pressure is higher than the value inferred under the assumption of 
thermal equilibrium. Markevitch et al. (1996) point out that this possibility might be real since 
the equilibration timescale is comparable to the typical time since the last merger. 

In principle, the deviation from thermal equilibrium could result either from smooth accretion 
or from mergers of clumps. For simplicity, we focus on smooth accretion in this paper. Making use 
of the versatility of a semi-analytic approach, we employ a self-similar model of spherical accretion 
to examine the regime of validity of the assumption of thermal equilibrium between the electrons 
and ions. This approach is similar to that of Itoh (1978), who considered thermal equilibration in 
supernova remnants using a self-similar blast-wave model. We assume that Coulomb collisions are 
the only coupling mechanism between the temperatures of electrons and ions. It can be shown 
that a homogeneous plasma with isotropic, monotonically decreasing distribution functions is 
stable against plasma instabilities (Rosenbluth 1965). Additional coupling might be provided by 
plasma instabilities in more complicated conditions, but we ignore them in this discussion. 

In §2 we describe the method of the calculation. In §3 we relate the parameters of our model 
to the conditions in typical clusters. In §4 we describe the details of the numerical calculations. 
In §5 we quantify the effects of the approximations made in our model. Finally, §6 discusses the 
results and the possibility of detecting the difference between the ion and electron temperatures 
with future X-ray telescopes. Throughout the paper we assume a density parameter H = 1, a 
Hubble constant /150 = 1, and a hydrogen mass fraction X = 76%. 


2. Method of Calculation 
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2.1. Collision Timescales 


We first consider collisional relaxation processes in plasmas consisting of several species 
of particles. Spitzer (1962) showed that each species achieves a Maxwellian distribution on a 
timescale of order the collision time, 

my^(3/cr)3/2 (r/10®K)3/2 ^1/2 

“ S.Tlvrne^ZMnA “ (n/10-3cm-3) ZMnA’ 

where A = ruxlnip is the particle mass of the species x in units of the proton mass, Ze is the 
particle charge, and In A is the Coulomb logarithm. For temperatures above 4 x 10® K, 

In A PS 37.8 + ln(T/10® K) — - ln(n/10“® cm“®). (2) 

Thus, the collision time for protons, tpp = (mp/me)^/^ tee; is longer by a factor of 43 than that for 
electrons, tee- 


Spitzer also found that the exchange of energy in collisions between test particles and held 
particles, with Maxwellian distributions of temperatures T and 7/, is governed by the equation 


(IT _Tf-T 
dt teq 


( 3 ) 


where the equilibration timescale, feq, depends on the related particle masses, m and m/, and 
charge numbers, Z and Zf, as well as the density of held particles, nj, and is given by 


t 


eq 


3mm f (kT 

8(27r)V2njZ2z2e4 In A \ m ^ ~r^) 

1 86 X ipSyr (T/10«K)3/^ 

^ Z2z2AV2lnA(n//10-3cm-3) 



( 4 ) 


For collisions between protons and ions of similar mass, such as He"*"^, the equilibration timescale 
is comparable to tpp. For electrons and ions, the equilibration timescale is greater than tpp by a 
factor of order {mplme)^^"^■ 


Thus, equilibrium in an electron-proton plasma is achieved in several stages. First, the 
electrons reach a Maxwellian distribution with temperature Tg on a timescale tee- Then, on 
a longer timescale, tpp ~ {mpjme)^^‘^tee, the protons reach a Maxwellian distribution with 
temperature Tp. Finally, the two temperatures equalize on a timescale of order tep ~ {mp/me)tee- 
We would like to consider a fully ionized plasma of hydrogen and helium, which is only slightly 
more complicated. All three ion timescales, t^^++p, tpp, and tHe++He++) much shorter than the 
electron-ion equilibration timescale. Therefore, when considering electron-ion equilibration, we 
can regard the ions as sharing a single temperature Tj. In this case, 


dTe _ Ti - Te 

dt tei 


( 5 ) 
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where the energy exchange rate between the ions and the electrons is just the sum of the rates due 
to exchange with the two ion species, 

1/tei = 1/tep + l/^eHe++- (6) 


We will generally be in the regime in which Te/Tj ^ me/vfip, where the contribution from 
each species to t~l is proportional to UfZj/Af. For both p and He"*"^, this factor simply equals 
the field particle density, Uf. Thus, in the case of a fully ionized plasma of hydrogen and helium. 


tei = 7.97 X 10® 


(re/10®K)3/2 

vr--— - --- 

(ni/10“^ cm“^) In A 


(7) 


where n* is the total ion number density. This timescale differs from tep only by a factor of Up/rii. 


2.2. Self-Similar Model for Spherical Accretion 

For simplicity, we consider spherical accretion of gas onto an X-ray cluster in an 17 = 1 
universe. Bertschinger (1985) found that an initial top-hat mass perturbation results in a 
self-similar solution, for both purely collisional and purely collisionless fluids, as well as for a single 
collisional fluid in a potential well dominated by a collisionless fluid (i.e. <C 17 = 1.) 

We would like to adapt this last model to the case of a fully ionized plasma of hydrogen and 
helium, under the gravitational influence of collisionless dark matter. The gas in Bertschinger’s 
model is considered adiabatic, and we retain this simplifying assumption. Thus, we neglect 
radiative cooling, a well-justified assumption for the outer parts of clusters where the free-free 
cooling time is much longer than the Hubble time. We also assume that there exists a magnetic 
field which is sufficiently strong to suppress thermal conduction, yet small enough to have a 
negligible dynamical effect. These conditions are indeed satisfied by the inferred magnetic field 
strength ~ 0.1 — 1 pG in cluster environments (Kim et al. 1990; Kim, Tribble, & Kronberg 1991; 
Rephaeli 1988). The general case of a stationary shock in a fully ionized plasma with conduction 
is complicated by the fact that the electron thermal speed exceeds the shock speed, allowing the 
electrons to preheat the plasma ahead of the shock (Shafranov 1957; Zel’dovich & Raizer 1967). 
If conduction is suppressed by a magnetic field, however, we will see that electron-ion equilibrium 
can be treated entirely locally, which simplifies the problem considerably. 

The self-similar model is divided into three regions; pre-shock, shock, and post-shock. In the 
pre-shock region, the gas is taken to be cold, so that pressure is neglected. In our model, we take 
the plasma to be fully ionized even before the shock, so it cannot literally have zero temperature 
(e.g., photoionization by the UV background would heat the gas to 10^ ®K). However, as 
long as the temperature is sufficiently low that the initial pressure is dynamically insignificant 
compared to gravity (i.e. T <C 10^ K for X-ray clusters), this will not affect the pre-shock behavior. 
In this case, the acceleration is due to gravity alone, and thus the pre-shock boundary conditions 
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for the gas density, pressure, and velocity {p, p, and v) are the same as in the case of a single, 
cold, neutral gas. 

The behavior of the gas at the shock is determined by the shock jump conditions for p, p, and 
V. These conditions derive from the conservation of total mass, momentum, and energy. Since 
we neglect conduction, the equations relating the pre-shock density, pressure, and velocity, pi, pi, 
and vi, to the post-shock values, p 2 , P 2 , and V 2 , are the same for an ionized plasma as for a single 
monoatomic gas. Since the pre-shock flow is cold, pi is negligible compared to p 2 , and the shock 
is strong. 

Even though the relation between pre-shock and post-shock conditions is unchanged, the 
self-similar solution could still be altered by the finite shock thickness. The thickness of collisional 
shocks is typically of order the mean free path, which, for Coulomb collisions, is proportional to 
T^. For the ions, whose post-shock temperature is high, the mean free path in the post-shock 
region may be large. In the hottest clusters, at z = 0, it can become a non-negligible fraction of the 
shock radius, r^, which characterizes the scale of the solution. This could break the self-similarity 
of the solution. In the presence of the magnetic field mentioned above, however, a collisionless 
shock, which can be much thinner than a mean free path, is more likely. We will neglect the finite 
thickness of the shock in our model. 

The shock jump conditions determine the mean temperature, T, given the mass per particle, 
pnip. However, in order to find Tj and separately, we must consider the shock region in more 
detail. Inside the shock, collisions convert bulk kinetic energy into thermal energy. The increase, 
|A:AT, in thermal energy per particle is of order m{vi — ^ 2 )^. Because of the ratio of masses, AT 
is clearly much greater for the ions than for the electrons. Adiabatic compression increases the 
electron temperature by a factor of {P 2 /Pi)^~^, where 7 is the ratio of specific heats. However, 
P2/P1 = (7 + l)/(7 ~ 1) for a strong shock, so for 7 = 5/3, the electron temperature increases 
only by a factor of 2.5. Thus, as long as the pre-shock temperature, Ti, is negligible compared to 
the post-shock temperature, T 2 , adiabatic compression will not produce a value of Te 2 which is 
significant compared to Tj 2 . There could also be collisional energy exchange within the shock. As 
long as the shock is thin, the time spent in the shock will be much less than tgiiTg = T). When 
Te <C T, tei is shorter by (Tg/T)^/^, but we will see in §5 that this has only a small effect at later 
times. We therefore assume that immediately after the shock, the ions have all of the thermal 
energy, and that Te 2 = 0 . 

In a thin shock, it is possible that the ion species will not have time to relax to Maxwellian 
velocity distributions or to equilibrium with each other through collisions. However, since the 
contribution per ion to is the same for p and He"*"^, dT^/dt depends only on the mean ion 
temperature. Furthermore, once the ions reach equilibrium with each other, the rate of change of 
their temperatures due to collisions with electrons is the same, so they will remain in equilibrium. 
Thus, for simplicity, we assume a single ion temperature, Tj. 


In the post-shock region, the adiabatic equation must be replaced by an energy equation for 
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each species, describing dT/dt due to collisions as well as adiabatic compression. Assuming a 
single ion temperature reduces the number of energy equations to two; (i) an electron equation, 

dTf. Ti — Tf. Tf. dn 

-IT = 7-1 

dt tei n dt 

describing the temperature evolution of electrons in a given fluid element, and (ii) an equation for 
the mean temperature of the particles in the same fluid element. 


dT , dn 

— = (^_ i)_ 
dt n dt 

which is obtained by averaging the equations for the individual species. 


(9) 


In principle, there are also separate continuity and momentum equations for each species. 
However, the Coulomb force prevents charge separation on scales larger than the Debye length, 
and spherical symmetry precludes electrical currents without charge separation, so the electron 
and ion velocities are identical, and separating these equations provides no additional information. 
Furthermore, the total pressure, p, which appears in the total momentum equation, depends only 
on the mean temperature, T, which is independent of the extent of electron-ion equilibration and 
is governed by the same adiabatic equation as in the case of a single collisional gas. Thus, the 
fluid equations for the density, pressure, and velocity are unchanged, and are merely supplemented 
by equation which describes the temperature equilibration. This is important, since the 
equilibration timescale would otherwise break the self-similarity of the solutions for p, p, and v. 

Since, in the zero conduction case, the equations and boundary conditions determining p, 
p, and V in all three regions of the self-similar model are unaffected by the replacement of a 
single collisional gas with a fully ionized plasma, we obtain the exact same solution for these 
variables. Thus, we can regard Bertschinger’s solution as a background in which we can solve for 
Te according to equation (^). 

The self-similar solution has a single length scale, which we take to be the turnaround radius, 
rta{t). In the case of a top-hat perturbation, rta oc t®/®. We adopt the dimensionless functions 
from Bertschinger (1985) and write p, p, and v as functions of A = rjrta- 

p{r,t) = pcVtbD{\), 

p{r,t) = -P(A), 

v{r,t) = ^H(A). (10) 

The scaling in terms of the critical density pc = (OvrCf^)”^ follows from dimensional analysis, 
since far outside the initial top-hat perturbation, there is no other density scale. The factors of 
are convenient in the Db 1 case. When we consider the collisionless dark matter generating the 
potential well, we will also need the total mass of the system 

47r Q 

mtot{r,t) = —pcrl^M{X). 


( 11 ) 



- 7 - 


2.3. Electron-Ion Eqnilibration 


It is useful to introduce scaled temperatures, Tg = Tg/T and T* = Ti/T. We then find 

dTg 1 dTf, Tg (IT Ti — Tg 
dt T dt dt tei ’ 


where the adiabatic compression terms cancel. Note that tei still depends on Tg = TTg. However, 
under the post-shock adiabatic conditions, T oc cx; , so by coincidence the factor jui 
is actually a constant for a given fluid element. There is still some residual dependence in the 
Coulomb term, A, but since tei depends only logarithmically on A, this dependence can be 
neglected. We consider any fluid element just after the time, when it passed through the shock. 
It is convenient to define t 2 s to be the value which tei would have had just after the shock if Tg 
were unity, i.e. 


= tei{ts)Te 


7.97 X 10^ yr (T2(ts)/10^ K)^/^ 

lnA2g (ni2(ts)/10“3cm-3) 


(13) 


where A 2 S is, similarly, the value A would have had just inside the shock radius, if Tg were 1. 
Then, neglecting the slight variation in In A, 


and 


Finally, by definition, 


tei = t2.rg3/2 

(14) 

dTe Ti — Tg 

dt ~ t2sfl'’^ ' 

(15) 

— _ UeTe + UiTi 
~ Ue + rii 

(16) 


and by eliminating Tj from equation (15), we find 


1 /n^\ 3/2 

dt t 2 s \ rii J 


(17) 


For (1 — Tg) <C 1, this simply represents exponential decay of the fractional temperature 
difference, Tj — Tg, with an e-folding time of t 2 sni/{ni + rig) ~ t2s/2. In fact, equation dO) can 
be integrated analytically for arbitrary Tg to find t(Tg). If we neglect the pre-shock electron 
temperature as well as energy exchange within the shock, and assume Tg 2 = 0 at the time tg when 
a particular fluid element has just passed through the shock, then 


At = t-ts 


t2s 


rij \ 
m + UeJ 



(18) 


This solution has precisely the same form found by Shafranov (1957) in the case of a stationary 
planar shock without post-shock adiabatic compression. A nonzero value of Tg 2 would be 













equivalent to adding a correction, 6t, to At. Because of the rapid relaxation when Tg <C 1, this 
correction is generally <C t 2 s- Thus at late times, when Tf, ^ Te 2 , the corresponding correction to 
Te{t) is small. We will estimate the magnitude of this effect in §5. 


2.4. Cluster Profile at a Given Time 


We now have an implicit solution for the scaled electron temperature, T^, of a fluid element as 
a function of the time since that fluid element passed through the shock. For application to X-ray 
clusters, one is more interested in Tf> as a function of radius at a fixed cosmic time, t. We consider 
a fluid element observed at radius r and ask at what time, tg, it passed through the shock. The 
shock is characterized by a fixed value. As, of A = r/rta- For a given fluid element. 


dX ^ r ^ r rta 
dint rta rtaTta 


(19) 


Expressing r = v In terms of the dimensionless V = vt/rta and noting that trta/fta = 8/9 for infall 
onto a top-hat perturbation, we can rewrite equation (H) as 


dX 

dint 


= V--X, 
9 ’ 


( 20 ) 


so that 

t dX , , 

The function e(A) is 1 for A = As and increases as A decreases. For a fixed t, fluid elements at 
different r, and hence different A, have different values of tg. We are actually interested in 


At = t — tg = tg{e — 1) = t{e — l)/e. 


( 22 ) 


The timescale for relaxation is fixed by t 2 g, which is a constant for each fluid element but 
varies from one fluid element to the next. Since t 2 g oc T 2 we see from the scaling of p 

and p at fixed A (eq. Q| ) that 


t2s oc 



oc oc e-5/3. 

t = ts 


(23) 


Following equation (^), the scaled temperature, Tg, is determined by the ratio 


At 

t2s 


- 1 ). 

12 


(24) 


where t 2 is the value of t 2 s for the fluid element passing through the shock at time t. Combining 
this with the implicit solution for Te(At), and equation (21), we have an implicit solution for 
Tg as a function of A < A*. Given ng/uj, Tg(A) is parameterized by the value of the ratio t/t 2 - 
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These solutions are shown in Figure 1 for various values of this parameter, using Ue/nj = 1.07 as 
appropriate for a fully ionized hydrogen-helium plasma with X = 76%. For t/t 2 S> 1, thermal 
equilibrium takes place rapidly compared to the dynamical time, and the temperature difference 
is significant only in a narrow region near rghock; where the fluid elements have passed through 
the shock very recently. When t/t 2 ^ 1, the effect may be significant over a larger fraction of 
the shock radius. The detailed shape of the curves depends on e(A), which is determined by the 
velocity of the fluid elements as they settle toward their equilibrium radii. 


3. Model Parameters 


To determine the significance of the lack of electron-ion equilibrium, we need to relate the 
model parameters to values for typical observed clusters. The self-similar model itself has only one 
free parameter, which determines its scale. The mass, density, and temperature all depend on 
rta- Unfortunately, X-ray measurements of these quantities tend to come from the inner parts of 
clusters, where the approximations leading to the self-similar solution break down and the solution 
has unrealistic divergences. We need to relate an observable quantity in the inner regions of a 
cluster to a quantity in the outer regions where self-similarity holds. One way to do this is to use 
a simplified model, such as that considered by Eke, Cole &: Frenk (1996), in which both the X-ray 
gas and total mass have the form of a singular isothermal sphere in the virialized section of the 
cluster. In this model, the mass within a mean overdensity of 178 and the temperature of the 
X-ray gas are related by 


mi78 = 2.34 X IO^^/iso^Mq 


\( ^ 

1 1 

AlO^Ky^ 

(1 + ^)J 


I 3/2 


(25) 


where z is the cosmological redshift, and where we have assumed X = 76%. While this model is 
obviously an oversimplification, the relationship it predicts between mass and emission-weighted 
X-ray temperature agrees well with N-body/hydrodynamic simulations (Navarro et al. 1995; Pen 
1997). Furthermore, it meets our need to relate an observable quantity to a property of the outer 
parts of clusters, where self-similarity might hold. 


It is easy to show that the overdensity in Bertschinger’s collisionless model, given by the 
dimensionless quantity M(A)/A^, equals 178 at A ~ 0.255 where M ~ 2.98. Using the scaling 
relation equation (pT|), we find 


'^ta 


13.2/15-0' 

14.0 


Mpc(l -h z)~^ y 

Mpc(l + ^)-3/2 


W,178 


2 X 1015 hXMa 


'T \ 1/2 
-^x \ 


108 kJ 


1/3 


(26) 
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The shock occurs at a fraction As = 0.3472 of the turnaround radius, so 


Ts = 4.60 Mpc {1 + z) 




mi78 


1/3 


^2 X 10^® , 

/ T \ 

= 4.84AS,‘Mpc(l + zr»''2(^) . 


(27) 


Given rta and t = to(l + we can convert the dimensionless quantities H, P, and V of 

the self-similar solution to the physical quantities p, p, and v. In the self-similar model, the 
temperature T on P/D, which diverges as at small A. Since we do not expect the self-similar 

model to be realistic at small radii, we simply cut it off at A = 0.105 where T becomes greater 
than Tx. For typical cluster conditions, the electrons and ions reach thermal equilibrium outside 
this radius, and our results depend only on the calibration of r^, not on the details of the central 
temperature profile. 


In order to calculate the equilibration timescale, t 2 , we need the mean temperature, T 2 , the 
ion density, ni 2 , and the Coulomb logarithm, A 2 , all evaluated just inside the shock radius. Given 
the mean molecular weight, pnip (which we take to he p = 0.59 for X = 76%), we can obtain T 2 
from P 2 IP 2 - Using Bertschinger’s values for P 2 and D 2 , we find 


Ta = 0.44rx. 


(28) 


Similarly, we can find ni 2 and ne 2 from p 2 . 


I ni2 

\ne2 



X 10 ® cm ^ 



(1 + 4)3. 


The Coulomb logarithm, which depends on T and Ue, is then 


In A 2 


39.7 In 


108 K 


- iln 


0.10 


|ln(l+z). 


Combining these values, we find 


t 

t2 


0.828 + 


rp \ -3/2 

iIk) 




(29) 


(30) 


(31) 


Measurements of deuterium abundance in quasar absorption systems indicate that 
Uf,/i 5 Q = 0.096lgQ2o (Tytler, Fan, & Buries 1996; Buries & Tytler 1997). We assume 
Ubhio = 0.10. 


4. Numerical Calculation 

In order to produce cluster profiles at a fixed time, we need to reproduce the self-similar 
solutions for D, P, and V. We integrated the fluid equations for one-dimensional self-similar 
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accretion in an external potential, following the same techniques as Bertschinger (1985). 
Specifically, we integrated Bertschinger’s equation (3.4a-c), with M replaced with the collisionless 
Mx appropriate for the flf, <C 1 case, using the change of variables given by Bertschinger’s equation 
(3.12) to avoid singularities at the origin. The initial post-shock conditions for the integration were 
given by Bertschinger’s equation (3.6), where Vi and Di, as well as the location, A* = 0.3472, of 
the shock, were taken from Bertschinger’s Table 10. For simplicity, we calculated the collisionless 
mass distribution, M^, by quadratic interpolation of Bertschinger’s Table 4, using the scaling 
M oc for A < 0.02. 

The equations were numerically integrated inward from the shock using a fifth-order 
Runge-Kutta routine. Values of D, P, and V were calculated on a grid with AA = A/500. The 
shock location was determined by the condition that the velocity, V, should be zero at the origin. 
Because we use values for interpolated from a table, we are effectively using a slightly different 
mass distribution for the collisionless fluid. Thus A^ = 0.3472 is no longer precisely the right value. 
This causes discrepancies with Bertschinger’s ri;, <C 1 solution which are worst near the origin. 
Since we cut off our solutions for the electron and ion temperatures at A = 0.105 where T = T^, 
only the region outside that radius concerns us. In this region, P/D, which determines T, differs 
from Bertschinger’s results by less than 0.5%. While the fractional errors in V reach 7%, the 
resulting fractional errors in e are only 0.5%. 

We also used the two local integrals of motion given in Bertschinger’s equations (2.27) and 
(3.10) as checks on the calculation. For A > 0.105, the largest fractional errors in these integrals 
were 10“^. 

In Figure 2, we show the electron, ion, and mean temperatures as a function of radius for 
various values of the emission-weighted X-ray temperature. In making comparisons with observed 
X-ray clusters, the density could in principle be measured directly, rather than calculated from 
equation (p^. 


5. Approximations 


In §2, we made a number of approximations which we now proceed to justify given the model 
parameters estimated in §3. Several of the approximations result in errors near the shock, where 
the temperature equilibration is rapid due to the small value of Tg. The effect of these errors at 
later times is equivalent to a change of definition of At from t — tg to t — tg — St. Because the 
relaxation rate for Tg ~ 1 is much smaller than when Tg <C 1, the error STg ~ 6t{dTe/dt) at later 
times is generally quite small. One such approximation was the assumption that Tg = 0 at 
neglecting the pre-shock electron temperature as well as energy exchange within the shock. In this 


case, St is the time it would take for Tg to go from zero to the actual Te 2 , which equation (18) 



- 12 - 


shows is 



Thus, the fractional error in Tg is 



Te 




(1 - re)r-5/2 


(32) 


(33) 


for Te 2 ^0.1. This correction is small except where Tg is of order Tg 2 - 

To determine the significance of this effect, we need to know Te 2 - For T 2 = 0.44Tx, adiabatic 
compression alone results in Tg 2 = 5.7Tgi/T^. Thus, in X-ray clusters with ^ 10^ K, the 
resulting Te 2 will be small unless the pre-shock temperature is ^ 10® K. 

The contribution of collisional energy exchange within the shock to Te 2 depends on the 
detailed shock structure, but we can estimate it. If we repeat the derivation of §2.3, but ignore 
adiabatic compression and scale the temperatures in terms of r 2 instead of T, the factor of 
(1 — Tg) in equation (p)7|) is replaced with T/T 2 — Tg. This reduces dTg/dt by a factor > T 2 /T. 
The resulting Tg 2 is given by equation (|3^, where now 5t is the time spent in the shock multiplied 
by the time-average of T/T 2 . As a rough estimate, suppose this effective time, 5t, equals the 
post-shock proton collision time, tpp 2 , which results in Tg 2 = 0.66. Figure 3 compares the resulting 
electron temperature profile with the St = 0 result for = 10® K. The fractional difference is large 
at the shock, but quickly becomes smaller than the difference between Tg and T. As a fraction 
of (1 — Tg), the effect of energy exchange within the shock is less dramatic for lower temperature 
clusters. We should note that for St = tpp 2 , the width of the shock is of order 1 Mpc x (Tx/10® K)^ 
at 2 ; = 0, which in the hottest clusters is a significant fraction of the shock radius. This will 
break self-similarity and alter the overall density and temperature profiles of the cluster. Figure 3 
should only be regarded as an approximate description for the electron temperature in this case. 
Conversely, as long as a thin shock is a reasonable approximation (e.g. due to embedded magnetic 
fields), collisional energy exchange within the shock will not significantly alter the results for the 
electron temperature over most of the region where Tg and Tj differ. 


A second approximation which may fail near the shock is the assumption that Tg/Tj S> nig/nip. 
For Tg ^ 10“®, this breaks down, so the equilibration timescale levels off, rather than continuing 
to decrease as Tg^'^, and Tg grows only exponentially, rather than as oc (At)^/®. If Tg 2 is less than 
this value, our assumption results in a shift, St, opposite in sign to the one above. Given Tg 2 
estimated before, this is likely to be unimportant. Even for Tg ^ 10“®, tgi ~ 10“'^t2s; so the time 
\St\ to go from Tg ~ 10® K to Eg = 10“® is less than 4 x 10“'^t2s- Thus, even in this worst case, STg 
remains significant only as long as Tg ^ 10“®. 
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A final approximation which causes errors near the shock is our use of Te = 1 in lnA 2 s. The 
error introduced depends on the value of lnA 2 s, formally breaking the one-parameter family of 
solutions shown in Figure 1. However, for the values of lnA 2 s calculated in §3 and shown in 
Figure 2, we can estimate its magnitude. A rigorous upper bound on this effect can be calculated, 
and for Te 2 > 1000 K, it produces 5Te/Te < 5% for Te > 0.1. 

Two other approximations also involve In A, but are exact at the shock and become less 
accurate at smaller radii. These effects also break the one-parameter family of solutions and 
affect both Figure 1 and Figure 2. First, our scaling for t 2 s for a particular fluid element, given 
in equation (^), ignores the variation of A with T and n as the fluid element is adiabatically 
compressed. Since A T/Vn r\j T^^^ for adiabatic compression, and since T changes only by a 
factor of 2.2 from the shock to the centrally-dominated T^, the fractional effect on In A is less than 


0.6% and the integrated effect is even smaller. Second, the scaling used in equations (^) and 
to relate t 2 s to t 2 ignores the variation in lnA 2 s between fluid elements, since it replaces lnA 2 s 
with lnA 2 . However, for the values of Tx used in Figure 2, the maximum fractional error in Tg is 
0 . 1 %. 


6. Discussion 


Our model predicts significant differences between the electron and ion temperatures in 
regions extending up to a third of the way inward from the accretion shock. The effect is greatest 
in the hottest clusters, where the equilibration time, tgi, is longest. Because the equilibration 
is exponential (except at early times) with a timescale neither much longer nor much shorter 
than the ages of clusters, it is difficult to predict the extent of nonequilibrium precisely. As more 
direct information about the physical characteristics of the X-ray gas in the outer parts of clusters 
becomes available, it will be possible to refine these predictions. 


In principle, the ion temperature could be inferred by measuring the width of X-ray lines 
broadened by the Doppler shift due to the ion thermal velocity dispersion. The equilibration 
timescale, tgq, is proportional to A/Z'^, so the high ionization states which dominate at cluster 
temperatures should equilibrate with the protons on timescales < tpp. Thus we can assume a 
single ion temperature, Tj. The rms velocity dispersion due to thermal motions along the line of 
sight is vios = {kTi/rrii)^^'^, which produces a line with an rms width 

_ vios _ 1 [Wj _ 3.0 X 10-3 / Tj 
E c c\ mi VA \W^k) ^ ' 


where A = mi/mp. In terms of a resolving power, E/AE, where AE is the Full Width at Half 
Maximum of a Maxwellian velocity distribution, equation (|^ corresponds to 


E 

'ae 


1.4 X IO'^VA 




1/2 


(35) 
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For the iron lines near 7keV, a resolving power of E/AE ~ 1000, or AE = 7 eV, would be 
needed. The microcalorimeters of the XRS instrument on the planned Astro E satellite have 
an instrumental AE ~ 11.5 eV. Convolution with the intrinsic width would broaden this to an 
observed AE of 13.5 eV. A high signal-to-noise ratio, as well as an accurate calibration of the 
instrumental profile, would be necessary to detect this broadening. 

The HETG spectrometer on AXAF will reach the necessary resolving power of 500 - 1000 
at energies ^ 2keV, for observations of point sources. It might therefore be able to measure the 
widths of absorption lines of O VIII, Ne X, or Mg XII in the spectra of quasars located behind 
the cluster. While the ionization state of these elements is determined by rather than Ti, 
the populations of these states are still likely to be small in the hotter clusters where we expect 
significant differences between Tg and Tj. 

Turbulent line broadening may also be significant in the outer parts of clusters. The ratio of 
thermal to turbulent broadening is 

'^los 
Vtmh 

where cs is the sound speed. Widths of several lines of species with A ^ {cs/vtuvhY would be 
necessary to distinguish thermal broadening from turbulent broadening. 

Fine structure lines in the ultraviolet might provide another probe of the ion temperature. 
At the temperatures in question, proton impact excitation can contribute significantly to the 
population of the upper levels of the Fe XVIII A974 and Fe XXI AI354 lines (Raymond 1997). Its 
effect on these lines could therefore be used to measure the ion temperature. 

Differences between electron and ion temperatures might be even more significant in 
clusters formed by mergers than in our simple, spherical accretion model. In a merger, the time 
since gas was shocked would be smaller, but preheating of the gas would be more important. 
Three-dimensional hydrodynamic simulations would be essential in exploring the extent of 
deviations from thermal equilibrium in this case. 
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Fig. 1.— Scaled electron temperature, Tg/T, as a function of scaled radius, r/rghock) for various 
values of the ratio between the cosmic time t = to/{I + z) and the electron-ion equilibration time 
at the shock, t 2 - 
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Fig. 2.— Electron temperature (solid), ion temperature (dotted), and mean temperature (dashed) 
as a function of radius for four clusters with emission-weighted X-ray temperatures, T^, of (a) 10® K, 
(b) 8 X lO’^K, (c) 6 X lO’^K, and (d) 4 x lO’^K. 









Temperature (K) 


- 18 - 


1X10 

8x10 

6x10 

4x10 

2x10 



1 2 3 4 5 6 

radius (Mpc) 


Fig. 3.— Effect of energy exchange within the shock for a cluster with emission-weighted 
temperature = 10® K. The electron temperature (solid) assuming an effective time within 
the shock of one post-shock proton-proton collision time {6t = tpp 2 ) is compared with the case of no 
exchange {St = 0). The mean temperature (dashed) is also shown. The effect is smaller for lower 
temperature clusters. 




